COMPUTING THE SUPPORT OF LOCAL COHOMOLOGY MODULES 



JOSEP ALVAREZ MONTANER AND ANTON LEYKIN 

Abstract. For a polynomial ring Ft = k[xi, ...,x n ], we present a method to compute the characteristic 
cycle of the localization Rr for any nonzero polynomial / 6 R that avoids a direct computation of Rt 
as a D-module. Based on this approach, we develop an algorithm for computing the characteristic cycle 
of the local cohomology modules H^(R) for any ideal ICR using the Cech complex. The algorithm, 
in particular, is useful for answering questions regarding vanishing of local cohomology modules and 
computing Lyubeznik numbers. These applications are illustrated by examples of computations using 
our implementation of the algorithm in Macaulay 2. 



1. Introduction 

Let k be a field of characteristic zero and R = k[x\, . . . , x n ] the ring of polynomials in n variables. For 
any ideal IC J?, the local cohomology modules Hj(R) have a natural finitely generated module structure 
over the Weyl algebra A n . Recently, there has been an effort made towards effective computation of these 
modules by using the theory of Grobner bases over rings of differential operators. Algorithms given by 
U. Walther and T. Oaku and N. Takayama [201 provide a utility for such computation and are both 
implemented in the package D-modules 16 for Macaulay 2 |10| . 

Walther's algorithm is based on the construction of the Cech complex in the category of A„-modules. 
So it is necessary to give a description of the localization Rf at a polynomial / e R. An algorithm 
to compute these modules was given by T. Oaku in |19) . The main ingredient of the algorithm is the 
computation of the Bernstein-Sato polynomial of / which turns out to be a major bottleneck due to its 
complexity. 

To give a presentation of Rf or Hj(R) as A„-modules is out of the scope of this work. Our aim 
is to provide an algorithm to compute an invariant that can be associated to a finitely generated A n - 
module, the characteristic cycle. This invariant gives a description of the support of the A n -module 
as an i?-module, so it is a useful tool to prove the vanishing of local cohomology modules. Moreover, 
the characteristic cycles of local cohomology modules also give us some extra information since their 
multiplicities are a set of numerical invariants of the quotient ring R/I (see Among these invariants 
we may find Lyubeznik numbers that were introduced in |17| . 

We will present an algorithm to compute the characteristic cycle of any local cohomology module. It 
comes naturally from the structure of the Cech complex and the additivity of the characteristic cycle with 
respect to short exact sequences. The requirement is that we have to compute first the characteristic cycle 
of the localizations appearing in the Cech complex. To do so, we present a method based on a geometric 
formula given by V. Ginsburg in 9 and reinterpreted by J. Briangon, P. Maisonobe and M. Merle in 0|. 
The advantage of this approach is that we will not have to compute the Bernstein-Sato polynomial of / 
and we will be operating in a commutative graded ring in 2n variables instead of operating in the Weyl 
algebra A n . The algorithm we will present is an elaboration of 4 Thm. 3.4.2], which - we have to point 
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out - is stated in the complex analytic context. For our computational purposes we are interested in the 
algebraic context since we will need (absolute) primary decomposition so we have to make sure that, at 
least for the examples we will develop, this result is also true in the algebraic counterpart. The complex 
algebraic case may already be found in [9] but the approach we will use in this work is through flat base 
change. It allows us to work over any field of characteristic zero for a large sample of examples since 
localization modules and local cohomology modules have a good behavior with respect to this operation. 
Since absolute primary decomposition is not implemented in Macaulay 2, we compute over the field of 
rational numbers. This replacement, however, is not an issue for the examples that we present. 

The scripts of the source codes we will use in this work as well as the output in full detail of the 
examples are available at the web page http://www2.math.uic.edu/^leykin/CC In the future, we would 
explore the possibility of using numerical primary decomposition - technique for varieties over C which 
is being developed by the second author. 

2. Basics on the theory of P-modules 

Let X — C™ be the complex analytic space with coordinate system x\, . . . ,x n . Given the convergent 
series ring R = C{x±, . . . , x n } consider the associated ring of differential operators D n := R(d\, . . . , d n ), 
i.e. the ring extension generated by the partial derivatives di = with the relations given by didj = 
djdi and dir — rdi = where r € R. For any unexplained terminology concerning the theory of rings 
of differential operators we shall use |3] , [7] . 

The ring D n has a natural increasing filtration given by the total order; the corresponding associated 
graded ring gr{D n ) is isomorphic to the polynomial ring R[a\, . . . , a n }. A finitely generated £> n -module 
M has an increasing sequence of finitely generated i?-submodules such that the associated graded module 
gr(M) is a finitely generated <?r(D n )-module. The characteristic ideal of M is the ideal in gr(D n ) = 
R[ai, . . . , a n ] given by the radical ideal J(M) := rad (Axm gr m n \(gr(M))). The ideal J(M) is independent 
of the good filtration on M. The characteristic variety of M is the closed algebraic set given by: 

C{M) := V(J(M)) C Spec {gr{D n )) = Spec (R[a x , . . . , aj). 

The characteristic variety describes the support of a finitely generated D„-module as i?-module. Let 
7r : Spec (R[ai, . . . , a n \) — > Spec (R) be the map defined by tt(x, a) = x. Then Supp H (M) = w(C(M)), 

We single out the important class of regular holonomic D„-modules. Namely, a finitely generated 
£> n -module M is holonomic if M = or dim C(M) = n. It is regular if there exists a good filtration on 
M such that Ann gr ( Dri }(gr(M)) is a radical ideal ([3], see also §3], [7]). 

The characteristic cycle of M is defined as: 

CC{M) = ^2rm A, 

where the sum is taken over all the irreducible components Aj = V(qi) of the characteristic variety C(M), 
where G Spec (gr(D n )) and m-i is the multiplicity of gr(M) at a generic point along each component 
Aj. These multiplicities can be computed via Hilbert functions (see H^dil)- Notice that the contraction 
of qi to R is a prime ideal so the variety 7r(Aj) is irreducible. These components can be described in 
terms of conormal bundles to Xi :— 7r(Aj) in A, i.e. 

CC(M) = ^m, T£.X. 

In particular, the support of M is Supp fl (M) = (JXj. For details we refer, among others, to |^ §10], 
HI §7.5]. 
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2.1. Characteristic cycle of a localization. Let M be a regular holonomic D„-module. Then the lo- 
calization Mf at a polynomial / £ R is a regular holonomic D„-module as well. A geometric formula that 
provides the characteristic cycle of Mf in terms of the characteristic cycle of M is given by V. Ginsburg 
in |5] and became known to us through the interpretation of J. Briancon, P. Maisonobe and M. Merle 
in El 

First we will recall how to compute the conormal bundle relative to /. Let Y° be the smooth part of 
a subvariety Y C X where /|y is a submersion. Set: 

W = {(x, a) £ T*X | x£Y° and a annihilates T x (/|y)- 1 (/(a;))}. 

The conormal bundle relative to /, denoted by Tt y , is then the closure of W in T*X\y. 

Theorem 2.1. (gl Thm. 3.4.2]) Let M be a regular holonomic D n -module with characteristic cycle 
CC{M) = m i T X .X and let f £ R be a polynomial. Then 

CC{M f )= J2 m ii V i+ T x i X ) 

with Ti = mijTij, where are the irreducible components of the divisor defined by f in T% x and 
rnij are the corresponding multiplicities. 

Remark 2.2. Assume for simplicity that M is a regular holonomic D n -module such that CC(M) = T Y X 
and let / £ R be a polynomial such that /(F) ^ 0. By the formula above we have CC(Mf) = TyX + F. 
It is worthwhile to point out that the reduced variety associated to T is the characteristic variety of the 
local cohomology module Hh^(M). 

Example 2.3. Set R = C{x, y, z}, M = HlJR), f = x and g = y. 

We have CC(R) = T* X X. Then Tt = {(x, y, z, a, b, c) £ T*X \ b = 0, c = 0} and the divisor defined 
by / in T*, is F = {(x, y, z, a, b, c) £ T*X \ b = 0, c = 0, x = 0} = Tf Q} X. Thus 

CC(R X )=T* X X + T* {X=0] X 

We have CC(M) = T* {x=Q} X. Then T* |{ _ q} = {(x, y, z, a, b, c) £ T*X \ c = 0, x = 0} and the divisor 
defined by g in T*^^ =g} is T = {(x, y, z, a, b, c) £ T*X \ c = 0, x = 0, y = 0} = T^ x=y=0} X. Thus 

CC(My) = T{ x=0 yX + T£ x=y=0 }X 

The multiplicities my appearing in the formula are the multiplicities of a generic point x along each 
component Tij of Tj and can be computed via Hilbert functions as in |15| . 

Lemma 2.4. Let e(F,x) denote the multiplicity of the variety T C T*X defined by the ideal I C 
R[a±, . . . , a n ] at a point x. Then, the multiplicity m of a generic point x along T is 

m = e(r,x)/e(T red ,x), 

where T red is the variety defined by rad (J) . 

Proof. A reformulation of ^2 Prop. 3.11] for the particular case of x being a generic point gives us the 
desired result, i.e. e(T,x) — e(T red ,x) ■ m. □ 

2.2. Algebraic P-modules. Let X = C" be the complex affiine space with coordinate system xi, . . . , x n . 
Given the polynomial ring R — C[xi, . . . , x n ] consider the associated ring of differential operators A n := 
R(di, . . . ,d n ), i.e. the Weyl algebra. The theories of algebraic P-modules and analytic I?- modules are 
very closely related. If one mimics the constructions given for the ring D ni one can check that the results 
we have considered before, conveniently reformulated, remain true for A n . In particular we may construct 
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an algebraic characteristic cycle as a counterpart to the analytic characteristic cycle described before. 
Our aim is to explain how both cycles are related. 

Set C{x} := C{xi, . . . , x n } and C[x] := C[xi, . . . , x n ]. Let M be a regular holonomic A„-module. 
The D„-module M an := C{x} ®c[x] M is also regular holonomic. For a good filtration {Mj},>o on 
M the filtration \Mf n := C{x} <8>c[x] ^i}i>o is a ls° good due to the fact that C{x} is flat over C[x]. 
Therefore gr(M an ) ~ C{x} ®c[x] 9 r (M) so the characteristic variety of M an is the extension of the 
characteristic variety of M, i.e. C(M an ) = C(M) an . However, we should notice that the components of 
the characteristic variety may differ depending on the ring we are considering. In particular we may have 
algebraically irreducible components that are analytically reducible. 

The regular holonomic A n -modules we will consider in this work, i.e. the polynomial ring R = C[x], 
the localization Rt for a polynomial / G R, and the local cohomology modules Hj(R), all have a good 
behavior with respect to flat base change. We state that, roughly speaking, the formulas of the algebraic 
and analytical characteristic variety of these modules are the same but the components and multiplicities 
of the corresponding characteristic cycle may differ. 

Remark 2.5. The results of this section can be stated in general for X being any smooth algebraic variety 
over C. It is worth to point out that M — > M an gives an equivalence between the category of regular 
holonomic T>x -modules and the category of regular holonomic 2?™-modules when X is projective (see 
13 §3]). 

3. Algorithmic approach to Brianqon-Maisonobe-Merle's Theorem 

From now on we will assume that R = C[xi, . . . , x n ] is the polynomial ring so we will be working 
in the algebraic context. Let M be a regular holonomic A n -module with algebraic characteristic cycle 
CC(M) —J2 m i T^.X and let / G R be a polynomial. Our aim is to compute the characteristic cycle of 
the localization M/ operating in the commutative graded ring gr(A n ) = R[a\, . . . , a n ]. We are going to 
provide two algorithms that are an elaboration of Theorem l2.ll The first one computes the part Ti of the 
formula in Theorem 12 . II corresponding to each irreducible component T^.X in the characteristic cycle of 
M. The second one computes the components and the corresponding multiplicities of the varieties Ti. 

Theorem 12. II is a geometric reformulation of a result given by V. Ginsburg Thm. 3.3]. Even though 
it is stated in the analytic context me may find in Ginsburg's paper the algebraic counterpart to the same 
result, see |H1 Thm. 3.2]. We may interpret it as in Section 2.2. through flat base change. 

Algorithm 3.1. (Divisor defined by f in the conormal relative to f) 

Input: Generators g%,...,gd of an ideal I C R defining the algebraic variety Y = V(I) C X and a 
polynomial / G R. 

Output: Divisor defined by / in the conormal Tt relative to /. 

Compute the smooth part Y° of Y where f\y is a submersion: 
(Oa) Compute V/ = ( J£) 

(Ob) Compute Y° = Y \ V(I°), where 1° C R is the defining ideal of {x G Y \ V/(x) = 0} and the 

singular locus of Y. 
Compute the conormal relative to / 
(la) Compute K — ker</>, where the <fi '■ R n ~ + R d+1 /I sends 

s^(Vf,Vgi,...,Vg d )-s€R d+1 /I. 

(lb) Let J C gr(A n ) — R[a 1: a n ] be the ideal generated by {(<ii, a n ) ■ b | b G K}. 
(lc) Compute J sat = J : {gr(A n )I°)°° . (Note: I{T* W ) = 

Compute the divisor defined by / in Ti, 
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(2a) Compute Kf = ker</>/, where the map <j>f : R n — ► R d+1 /(I + (/)) sends 

a -» (V/, Vfli, V 5d ) ■ « G i? d+1 /(i+ (/))• 
(2b) Let J/ C gr(A„) = i?[ai, a„] be the ideal generated by {(ai, a„) ■ 6 | 6 G if/}. 
(2c) C - J sat + (/) + J/ C .gr(A0- 
Return: The ideal C that defines the divisor f in T% 

Proof. (Correctness of the algorithm) The steps (Oa), (Ob) follow from the definition of f\y being a 

submersion. The relative conormal Tt, is the closure of 

j\y 

W = {{x, a) G T*X | x G F°,Vs G if, a(s) = 0}. 

For every point x EY°, the tangent space T X Y° is a specialization of V(K ), where K is computed in step 
(la). A defining ideal of W is produced in (lb) and, finally, taking the closure amounts to the saturation 
in (lc). In order to restrict to / = 0, it is not enough to compute J sa t + (/). However, step (2a) and (2b) 
that follow closely the idea of (la) and (lb) provide the necessary correction term in (2c). 

Recall that the analytic extension of the ideal C we obtain with the algorithm is what we would obtain 
applying Theorem 12.11 in order to compute the analytic characteristic cycle of the localization module 
(see Section 2.2). In our case, the ideal C will give us the components of the algebraic characteristic 
variety. 

□ 

Algorithm 3.2. (Components and multiplicities of the characteristic cycle) 

Input: The characteristic cycle CC(M) = J2 m i T^.X of a regular holonomic A„-module M and a 
polynomial / G R. 

Output: The characteristic cycle CC(Mj) = J2f(x i ) 7 to m i(Xi + T Xi X )- 

For every component Y — Xi we have to compute the ideal Cj corresponding to the divisor defined 
by / in T% using Algorithm 13. II Then: 

Compute the components of C, 

(la) Compute the associated primes dj of C\. 

(lb) Compute iy = dj (1 R (if you need to know the defining ideal of Xij — 7r(r,j) in Theorem 12. 1(1 . 
Compute the multiplicities 

(2) Compute the multiplicity in Theorem 12 .11 as the multiplicity of a generic point x along each 
component Cy of Cj as in Lemma 12.41 via Hilbert functions. 
Return: The components of CC(Mf) and their corresponding multiplicities. 

Proof. The correctness of the algorithm is straightforward and follows from Lemma l2~4l □ 

The algorithm we propose requires the computation of the associated primes of an ideal; primary 
decomposition is also needed in the implementation if we want to avoid choosing generic points when 
computing the multiplicities (see Lemma 12.411 . Therefore, we have to restrict ourselves to computations 
in the polynomial ring R = Q[xi, . . . , x n ] as we implemented the algorithm in the computer system 
Macaulay 2. What we are going to construct is the characteristic cycle of a regular holonomic A n - 
module where now A n stands for the Weyl algebra with rational coefficients. By flat base change we can 
extend the ideal C we obtain with Algorithm 13 . 1 1 to any ring of polynomials over a field of characteristic 
zero or to the convergent series ring over C. As we stated in Section 2.2, the primary components may 
differ depending on the ring we are considering. 

In order to construct the algebraic characteristic cycle over Q we would need to find the absolute 
primary decomposition of the ideal C we obtain with Algorithm 13.11 Even though the Macaulay 2 
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command for primary decomposition is not implemented over the algebraic closure of Q, it suffices for 
the examples we treat in the next section. 

Another fine point in the implementation is the treatment of embedded components of the ideal C 
outputted in Algorithm 13. II The ideal C contains the complete information about maximal components 
of the divisor / on Tf^, in particular, we can compute their multiplicities. However, the primary ideals 
in the decomposition of C that correspond to an embedded component may not lead to the correct 
multiplicity due to the global nature of our computations. In order to obtain this multiplicity we restrict 
the divisor to the embedded component, which amounts to rerunning Algorithm ^. II 'modulo' its defining 
ideal. The top-level routine in our implementation processes components recursively "descending" to, 
i.e., localizing at, the embedded components when needed. 

4. Characteristic cycle and Cech complex 

Let I — {fx, ... , f s ) C R = C[xi, . . . , x„] be an ideal and M be a holonomic ^4 ra -module. In this section 
we are going to compute the characteristic cycle of the local cohomology modules Hj(M) using the Cech 
complex 

s 

C'(fx, ...,f s ;M): • M - • M fi ^ ► M fl ... fs — ► 0. 

i=l 

For simplicity we will assume from now on that M is indecomposable. Otherwise, if M = Mi M2, 
then 

C'(h, ..../,: M) = C'(h, ..../,: Mr) 8 C'(fr, . . . , /,; M 2 ) 

and Hj(M) — Hj(Mi) iJJ(M 2 ) for all r, so we can compute the characteristic cycle of both local 
cohomology modules separately. Sometimes we will denote the localization modules appearing in the 

Cech complex Mf a , where f a — Yl a -=i /« f° r au a e (0) ^} s - We will also denote |ot| = a.\ H \-a 8 and 

£i, . . . , £ s will be the natural basis of If . 

^From the characteristic cycle of the localization modules in the complex we develop an algorithm 
to extract the precise information needed to describe the characteristic cycles of the local cohomology 
modules. The algorithm comes naturally from the structure of the Cech complex and the additivity of 
the characteristic cycle with respect to short exact sequences. However the following assumption will be 
required: 

(f) For all a £ {0, 1} S such that on — 0, the localization map Mf a — >Mf a+e . is either a natural 
inclusion, i.e., Mf a is saturated with respect to f Si , or Mf a+e , = 0. 

For unexplained terminology on the theory of complexes we refer to |22j. To shed some light on the 
process we first present the case of / being generated by one and two elements. 

4.1. The case s = 1. We have the short exact sequence 

— > H°j{M) — > M M h — ► H}{M) — >Q 

Under the assumption (f) either CC(Hj(M)) = CC(M h )-CC{M) if M is /i-saturated or CC(H^(M)) = 
CC(M) if Mf 1 = 0. The algorithm we propose to compute the characteristic cycle works in both cases 
and boils down to the following step: 

Prune the characteristic cycles of Mf t and M , where 'prune' means remove the components (counting 
multiplicities) that both modules have in common. 

The characteristic cycle of Hj(R) (resp. Hj(R)) is the formal sum of components of M (resp. Mf x ) 
that survived this process. 
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4.2. The case s = 2. The vertical sequences of the following diagram are exact 
C'(h;M): >M >M h >0 



C' (/i , / 2 ; M) : M M A M h M hh 

C* (/i ; M h )[-!]: M h M hh 

so we have an exact sequence of Cech complexes 

(i) — C*(/ i; M /2 )[-l] — > C"(/i, / 2 ;M) — <3»(/ i; M)— 

where [—1] stands for the result of shifting the complex one place to the right. Analogously, i.e. switching 
the roles played by f± and / 2 , we also have 

(ii) — C*(/ 2 — / 2 ;M) — C"(/ 2 ;M)^0 

Notice that the vanishing of any localization module reduces the Cech complex to the case s = 1. So 
we are going to consider the only case remaining under the assumption (f), i.e. M is saturated with 
respect to /i/ 2 . Consider the long exact sequence of cohomology modules associated to (i) 

0^H-^(M f2 ) — H°(M) H° fi) (M) H° fi) (M h ) Hj(M) ■ ■ ■ 
where & are the connecting maps. Non-vanishing may occur only in the sequence 

H}(M) H\ h) {M) H\ h) {M h ) H*(M) 
which breaks down into two short exact sequences with C\ = CokerJ 1 : 

— >H}(M) — >H^ fi) (M) — >d — >0 and — >d — >H^ fi) (M h ) — > Hj(M) — >0 

In order to get CC(H)(M)) and CC(Hj(M)) we only have to compute the characteristic cycle of C\ 
since we already know that CC(H^ h) (M)) = CC(M fl )- CC(M) and CC(H l (fi) (M h )) = CC(M flh ) - 
CC(M f2 ). 

Claim: CC(C\) — ^rrii T^.X where the sum is taken over the components (counting multiplicities) 
that CC(Hf fi) (M)) and CC(Hf fi) (M f2 )) have in common. 

Proof. Assume that there is a component T^.X in CC(HhJM)) and CC(HhJMf 2 )) not appearing 
in CC(d), i.e. T^X is a component of CC(Hj(M)) and CC{Hj{M)). This component shows up in 
the computation of the characteristic cycle of the cohomology of the subcomplex C"(/ 2 ; M/Jf— 1], since 
it is in fact a component of CC(Mf 1 ) and CC(Mf x f 2 ), but is not a component of CC(M) and CC(Mf 2 ). 
Consider the long exact sequence of cohomology modules associated to the complex (ii): 

>Hj(M) H° {h) (M) H° {h) (M fl ) Hj(M) H\ h) (M) ■ ■ ■ 

It follows that the component T^.X should belong to CC(H9 f AM h )) and CC(H} f s(M fl )) in order to 
fulfill the hypothesis of being a component of CC(H}(M)) and CC(Hj(M)). Thus we get a contradiction 
since iJ^^M/J = as Mf 1 has no / 2 -torsion. □ 

To summarize the computation of the characteristic cycle of Hj(M) using the sequence of Cech com- 
plexes (i) we propose the following algorithm: 

(1) Prune the characteristic cycle of Mf 1 and M. 
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Prune the characteristic cycle of Mf t f 2 and Mf 2 . 

(2) Prune the characteristic cycle of Mf 2 and M. 
Prune the characteristic cycle of Mf 1 f 2 and Mf ± . 

'Prune' means remove the components (counting multiplicities) such that both characteristic cycles 
still have in common in that step of the algorithm. The characteristic cycle of Hj(M) (resp. Hj(M), 
Hj(M)) is the formal sum of comand CC(Hhs (M/ 1 ))ponents of M (resp. Mf 1 and M/ 2 , M/^), i.e., 

components of th characteristic cycles of C°(fi, fa; M) (resp. C x {fx, fi\ M), C 2 (fi, fi\ M)) that survived 
to the process. 

Remark 4.1. Naturally, one may permute the steps (1) and (2) in the above algorithm. 

4.3. The general case. Let / = {fi, ■ ■ ■ , f s ) ^= R be an ideal and M be a holonomic A„-module 
satisfying 

(f) For all a <E {0, 1} S such that on — 0, the localization map Mf a — >Mf a+c . is a natural inclusion, 
i.e. Mf a is saturated with respect to f ei , or Mf a+E , = 0. 

Our aim is to proceed inductively in order to compute the characteristic cycle of the local cohomology 
modules Hj{M). To this purpose it is useful to visualize the Cech complex C'(fx, . . . , f s ; M) as a s- 
hypercube where the edges are the localization maps that, with the corresponding sign, describe the 
differentials of the complex. The exact sequence of complexes 







f "•;/:. . . . , / s _i; M f ,)[-1] C'(h, ..../,: .V/) 



can be easily identified in the s-hypercube. For the case s = 
C'{fu K h\ AO, C'ifu h\ M) and 6*(h,h\ M A )[-1] as follows 

M /l/2/3 



/2/3 




<?•(/!,... ,/ a _ i; M)—»0 
3 we visualize the Cech complexes 

M /l/2/3 




Chasing the diagrams, one may check that the localization map with respect to / s , i.e., the edges 
Mf a — *Mf a+e in the s-hypercube, induces the connecting maps S 3 in the long exact sequence of coho- 
mology modules 

^ H U u ...j^)( M fs) — H°(M) — > H° fi _ fs i) (M) H° fi _ fs i) (M f3 ) — H\(M) — > • • • 

We are not going to give a precise description of the connecting maps, since we are only interested in the 
data given by the characteristic cycle. The formula we obtain in Theorem 14. 41 for the characteristic cycle 
of the local cohomology modules Hj(M) is given just in terms of the components of the characteristic 
cycle of the localizations Mf a . The precise information we need to extract is given by the following 
algorithmic procedure. 

Algorithm 4.2. (Characteristic cycle and Cech complex ) 

Input: Characteristic cycles CC(Mf a ) = J2 m a,i T\.X for a € {0, 1} S . 
Output: (Pruned) characteristic subcycles CC(Mj a ) C CC(Mf a ) for a g {0, 1} S . 
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Prune the extra components 

For j from 1 to s, incrementing by 1 

(j) Prune the localizations Mf a and Mf a+e . for all a £ {0, 1} S such that ctj — 0, where 'prune' means 
remove the components (counting multiplicities) such that both modules still have in common 
after step (j — 1). 

Return: The components and the corresponding multiplicities of a characteristic subcycle of CC(Mf a ) 
for a G {0, 1} S . 

For the case s = 3 we visualize the steps of the algorithm as follows 

M fl f2 f3 M fl f2 f3 M /l f2 f3 




The solid arrows indicate the modules we must prune at each step. 

Remark 4.3. As in the case s = 2, the order we propose in the algorithm depends on the Cech subcom- 
plexes we consider when computing the characteristic cycle of the cohomology of the Cech complex. We 
can obtain equivalent algorithms permuting the generators of the ideal /. It is also worth mentioning 
that the algorithm can be used in the algebraic context over any field of characteristic zero and in the 
analytic context. 

Theorem 4.4. Let I = (/i, . . . , f s ) C R be an ideal and M be an indecomposable holonomic A n -module 
satisfying (f). Then 

CC(H r j(M)) = CC(M fa ), 

\a\—r 

where the pruned characteristic subcycles are obtained with Algorithm \4-S\ 

Proof. We proceed by induction on the number of generators s of the ideal. The cases s = 1, 2 have been 
already done. For s > 2 we have an exact sequence of complexes 

O — tf'CA,..., / s _i; M/,)[-l] — > C'(h, . . . , f s ; M) — ► C'{h, . . . , / a _ i; M)— »0 

Splitting the corresponding associated long exact sequence of cohomology modules 

^ H lfl..,f^)( M f.) — H "( M ) — ^(/ 1 ,...,/ 8 _ 1 )(^/J — W — ' ' ' 

into short exact sequences we obtain 

— >A r — >H r j(M) — >B r — >0 
0^B r -^H r {fi Jsi) (M)^C r ^0 
0— C r — »fq- /iiiii)/# _ i) (Af / J— »i4 r+ i— 

The characteristic cycle of i?L ^ (resp. i?L ^ ^(My s )) is the formal sum of components 

of Mf a satisfying a s = and |a| = r (resp. a s — 1 and |a| = r + 1) that survived to step (s — 1) of 
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the algorithm. Thus, for every r, in order to get CC(Hj(M)) we only have to compute CC(C r ) and use 
additivity of the characteristic cycle with respect to short exact sequences. 

Claim: CC(C r ) —J2 m i T^-.X where the sum is taken over the components (counting multiplicities) 
that CC(HL js(M)) and CC(HTf ^.(MfJ) have in common. 

Assume that there is a component TJ.X in CC{H r Uu ^{M)) and CC{H r (h Js i) (M fs )) not 
appearing in CC(C r ), i.e. it has not been pruned in step (s). In particular, this component belongs to 
the characteristic cycle of some localization modules Mf a and M/ Q+es satisfying a s = and \a\ = r. 
Then, this component is not pruned in the computation of the characteristic cycle of the cohomology of 
a convenient proper Cech subcomplex of C'{f\, . . . , f s ; M) containing Mf a and Mf a+Es . Thus we get a 
contradiction. 

□ 

If AI is not indecomposable we only have to apply Theorem 14.41 to each component. 

Example 4.5. Consider R — C[x] and the holonomic Ai-module M = R® H^(R). We have: 

CC(M) = T^X + T* {x=a] X 
CC{M x ) = T£X + T( x=0} X. 

Applying the pruning algorithm to each component, that satisfy (f), we get 

CC(H? a) (M)) = CC{Hl x) {M)) = T( X=0] X. 

One may be tempted to apply the pruning algorithm to M , however, it misleads us resulting in the 
seeming vanishing of the local cohomology modules H^(M). Notice that M is not saturated with respect 
to / = x. 

The question whether Theorem 14.41 still holds for indecomposable A n -modules not satisfying (f) is 
open. For some examples we may give an affirmative answer. 

Example 4.6. Set R — C[a;, y]. The holonomic ^-module M = H^ xy ^(R) is not saturated with respect 
to / = y. We have: 

CC(M) = T* x=Q} X + T* {y=Q} X + T* {x=y=0} X 
CC(M y )=T{ x=0} X + Tl x=y=0} X. 

Pruning the components they have in common we get CC(H^(M)) = T?_ Q yX. It agrees with the fact 
that H° y) (M)^Hf y) (R). 

Remark 4.7. Let / = (fi, . . . , f s ) C R be an ideal and M be an indecomposable holonomic A n -modulc 
saturated with respect to f± ■ ■ ■ f s . The first step of Algorithm 14 . 21 also comes from the fact that the Cech 
complex C*(/i, . . . , / s ; M) is quasi-isomorphic to 

'» • <> '' • M h /M M hf jM h ► M h ...fjM h ... fs — > 0. 

It would be interesting to continue the pruning process for the whole complex (not just for the com- 
ponents of the characteristic cycle) in order to find a complete description of a minimal complex quasi- 
isomorphic to the Cech complex. 

A canonical Cech complex is introduced in ^H] for the case of / being a monomial ideal and M = R. 
This complex is associated to a minimal free resolution of R/I in the same way as the usual Cech complex 
is associated to the Taylor resolution of R/I. The difference with the pruned Cech complex we propose 
is that, roughly speaking, they only prune when the localization map Rf a — >Rf a+e . is the identity. 
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5. Examples 

In this section we want to present some examples where we will apply our algorithm to compute 
the characteristic cycle of local cohomology modules. First we have to study localizations Rf of the 
polynomial ring R — Q[x\, . . . ,x n ] at a polynomial f € R. To compute its characteristic cycle directly 
one needs to: 

• Construct a presentation of the A^-module Rf, 
■ Compute the characteristic ideal J(i?/), 

• Compute the primary decomposition of J(Rf) and its corresponding multiplicities. 

The first two steps require expensive computations in the Weyl algebra A n since we have to compute 
the Bernstein-Sato polynomial of /. For some short examples we can do the job just using the Macaulay2 
commands Dlocalize and char I deal. 

Following the approach of this work we have developed some scripts written in Macaulay 2 that 
compute and print out the list of components and the corresponding multiplicities showing up in the 
characteristic cycles of the localizations Rf in the examples we present in this section. In fact we develop 
two different strategies that we may use depending on the examples we want to treat. 

• Single localization: Since the characteristic cycle of R is CC(R) = T^X, the characteristic cycle of 
Rf is CC(Rf) = T^X + r, where T is computed according to Theorem l2.1l so we may compute it in one 
step. Notice that the defining ideal of T may be quite large so computing its primary decomposition can 
be expensive. 

■ Iterative localization: We can apply Theorem 12.11 in an iterative way on the components of the 
polynomial /. This strategy is useful to treat large examples, since, usually it leads to computing 
primary decompositions of ideals of lower degrees compared to the former strategy. 

For the examples we present in this work both strategies can be applied. 

5.1. Local cohomology modules. Consider the ideal I C R = Q[xi, x@] generated by the minors of 
the matrix 

( xi x 2 x 3 \ 
\ x 4 x 5 x 6 ) ' 

It is a nontrivial problem to show that the local cohomology module Hj(R) is nonzero (see [121 
Remark 3.13], |13|). For example, Macaulay 2 runs out of memory before computing this module with 
the command localCohom. U. Walther |23l Example 6.1] gives a complete description of this module 
using a tailor-made implementation of his algorithm which is based on the construction of the Cech 
complex. The difference with the implementation of the Macaulay 2 command is that he uses iterative 
localization to reduce the complexity in the computation of Bernstein-Sato polynomials. 

Our method makes it possible to prove algorithmically that Hj(R) ^ from the computation of the 
characteristic cycles of the localization modules in the Cech complex which for this particular example 
looks like 

(*) -> R Rf, © R h © R h R hh © R hh © R hh -> R hhh 0, 
where fx = X\x§ — x 2 x^ f% — xiXq — X3X4 and fa = x 2 xq — X3X5. 

Remark 5.1. By flat base change we can also deduce the non- vanishing of the local cohomology module 
Hf(R) where R = k[x%, ■■■,xq] is the polynomial ring over any field k of characteristic zero. 

The list of components and their corresponding multiplicities showing up in the characteristic cycles 
of the chains in the Cech complex (★) and different from the whole space X that we get with our script 
contains 14 elements. A sample entry is as follows: 
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Component = V( ideal (xx -xx.xx -xx.xx - x x )) 

35 26 34 16 24 15 
entries-> HashTable{{0, 1, 2} => 2} 

{0, 1> => 1 

{0, 2} => 1 

-CO} => 

{1, 2} => 1 

{1} => 

{2} => 

Namely, the component corresponding to the ideal I is present with multiplicity one in Rf t f 21 Rf 2 f 3 
Rf 1 f 3 and with multiplicity two in Rf 1 f 2 f 3 - The following is the complete list of 14 components: 

M = V(fi), A 2 = V(h), A 3 = V(f 3 ), 

B 1 =V(x 3 ,x & ) 1 B 2 = V(x 2 ,x 5 ), B 3 = V(x 1 ,x i ), 

Ci = V(x 3 ,x 6 ,fi), C 2 = V(x 2 ,x 5 ,f 2 ), C 3 = V(x!,x 4 ,f 3 ), 

Di =V(x 1 ,x 2 ,x i ,x 5 ), D 2 = V(xi,x 3 ,X4,x 6 ), D 3 = V(x 2l x 3 , x 5 , x 6 ), 

E = V(x 1 ,x 2 ,x 3 ,x i ,x 5 ,x 6 ), F = V(I). 

Piecing the results of our computation together we can draw the 3-hypercube in Figure 1 . 




Figure 1 . Components of characteristic cycles for the Cech complex (*) (multiplicity > 
1 is specified in square brackets). 

To compute the characteristic cycle of the cohomology modules we have to apply Theorem 14.41 that 
has been implemented in the routine PruneCechComplexCC. According to the output 

O => {} 
{0} => O 
m => {> 
{2} => {} 

{0, 1} => -[ideal (xx - xx.xx - xx.xx - x x ) => 1} 
35 26 34 16 24 15 

{0, 2} => {} 
{1, 2} => {} 

{0, 1, 2} => {ideal (x , x , x , x , x , x ) => 1} 
6 5 4 3 2 1 

we get CC(H](R)) = T^X and C'C{H](R)) = T%X. Finally, it is worth to point out that the obtained 
result is coherent with the fact that the local cohomology module Hf(R) is isomorphic to the injective 
hull of the residue field Er(R/(xi, ...,xq)). 
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5.2. Lyubeznik numbers. Let R = k[xi, ...,x n ] be the polynomial ring over a field k of characteristic 
zero. Let / C R be an ideal and m = (xi, ...,x n ) be the homogeneous maximal ideal. G. Lyubeznik |17j 
has defined a new set of numerical invariants of the quotient ring R/I by means of the Bass numbers 

X p ,i(R/I) := fi p (m : H^-\R)) := dim fe Ext p R (k, H^\R)). 

These invariants can be described as the multiplicities of the characteristic cycle of the local cohomology 
modules Hl(H^ l {R)) (see 0). Namely, 

CC{Hl{Hr\R))) = X P ,i T* E X 

Lyubeznik numbers carry interesting topological information of the quotient ring Rj I as it is pointed 
in and [B]. To compute them for a given ideal ICR and arbitrary i,p we refer to U. Walther's 
algorithm 23, Algorithm 5.3] even though it has not been implemented yet. When / is a squarefree 
monomial ideal, a description of these invariants is given in pQ . Some other particular computations may 
also be found in |S] and |24j . 

Let I a R = Q[xi, ■■■,xq] be the ideal generated by the minors of the matrix 

Xi x 2 x 3 
X4 X5 Xq 

considered above, i.e. I = (X1X5 — X2X4,xiXg — X3X4,X2Xg — X3X5). We want to compute the characteristic 
cycle of the local cohomology modules Hm(H}(R)) for i = 2, 3 and Vp so we have to construct the Cech 
complex 

6 

(**) ^ M -h. M Xi ► M X1 ... X6 -> 0, 

i=l 

where M is either H](R) or Hf(R). Then we have to compute the characteristic cycles of the localization 
modules and use Theorem 14. 41 

• For M — Hf(R) we know that its characteristic cycle is T^X so, applying Theorem 12. II the Cech 
complex (**) reduces to the first term. Then, 

CC{Hl{Hj{R)))=T* E X 

and the other local cohomology modules vanish. 

• For M = H](R) we obtain 

CC(H 2 m (H](R)))=T E X 
CC(H 4 m (Hj(R)))=T E X 

and the other local cohomology modules vanish. We are not going to present the complete output with 
all the components as in Figure 1 for this case but at least we are going to show the multiplicities of the 
component T E X appearing in the Cech complex (**) in Figure 2. We point out that the components 
T E X that survive to Algorithm IO belong to CC(M XlX2 ) and CC(M XlXiX5X6 ). 

Using the properties that Lyubeznik numbers satisfy (see |171 Section 4]), we can collect the multi- 
plicities in a triangular matrix as follows: 

/0 1 0\ 

A(i?/J) = 1 



V V 
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M : M X1 


M X1X2 


E 




B[2] 


-^WxiX2X3X4 


£[3] 


^Mx 1X2X3X4X5 


£[3] M, 


M X2 


M X1X3 


E 


J\^x 1 x 2 x 4 


: £ 




£[3] 


^-^x 1X2^3^4^6 


£[3] 


M X3 


M xlXi 


: 


-^^xi X2X5 


: £ 


^^-x\X2X^xq 


E[3] 


^^xiX2X^x^xq 


£[3] 


M X4 


M XlXB 


E 


M^x\X2Xq ' 


£[3] 


MxiX2X4X& 


: E 


^^xiX2X^xqXq 


£[3] 


M X5 


M I1>6 


E 


-^^x 1X3X4 


: E 


■MxiX2X^xq 


E[3] 


^^X 1X3X4X5X5 


£[3] 


M X6 





E 


-^^x 1X3X5 ■ 


E[3] 


M-xix 2 X5Xq 


E[3] 


^$x 2 x 3 x 4 x 5 x e 


E[3) 




-^^X2^4 


E 


Af X ix 3 x 6 


: E 


M-Xl X3X4X5 


E[3] 








^x 2 x 5 


:0 


^^X 1X4X5 


: E 


■M-X\ X 3X4X5 


: E 








M X2 xq 




-^^X 1X4X3 


: E 


-^^X 1X3X5X5 


E[3] 












^xix 5 x 6 : 


E[3] 


-^xiX4X5Xg 


E{3] 








M X3X5 






E[3] 


^x 2 x 3 x 4 X5 


E[3] 








M X3X6 


:0 


JWx 2 x 3 X5 


: E 


^fx 2 X3X 4 x 6 


E{3] 








M X4X5 




^fx 2 x 3 x 6 


: E 


^fx 2 x 3 X5X 6 


: E 








M X4X q 




-^X2X4X5 


: E 


-^^X2X4X5Xg 


E[3] 








M X5X g 




-^^x 2 ^ 4 2; 6 ' 


E[3] 


3X4 £5 Eg 


E[3] 







^f X 2X§ X Q • E 

M X3XAX5 : E[3] 

^f X -$X4 X Q • E 
^^ X ^ X ^ X Q • E 

M X4 ,x 5Xe : E[2] 



Figure 2. Component T^X appearing in the Cech complex (**) (multiplicity > 1 is 
specified in square brackets). 

The complex variety V defined by / has an isolated singularity at the origin. The singular cohomology 
groups of V with complex coefficients and support at the origin can be described from Lyubeznik numbers 
(see In our case we get 

l = A 4 ,4 = dim c Hf 0} (V,C), 
1 = X 2 A=dim c Hf 0} (V,C), 
l = A , 3 =dim c tf { 3 0} (y,C). 

6. Conclusion and possible developments 

We have shown that characteristic cycles of local cohomology modules can be be computed by algorithm 
operating in commutative polynomial rings as opposed to the direct computation of these modules, which 
requires Grobner bases technique in noncommutative Weyl algebras. 

The computational engine of our method is primary decomposition, from which we extract only geo- 
metrical information. This prompts a natural interest in numerical primary decomposition, an algorithm 
that would produce just that - the descriptions of reduced components and their multiplicities - by means 
of numerical computations. 
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